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Abstract 

"With  air  traffic  movements  at  a  high  level,  techniques  to  assist  air  traffic  man¬ 
agement  using  computers  are  being  investigated.  One  technique  in  particular  being 
studied  is  the  early  adjustment  of  the  speed  of  arriving  aircraft  so  that  the  rate  of  flow 
near  to  the  airports  is  closely  matched  to  landing  capacity. 

A  Speed  Control  Adviser  has  been  developed  which  allocates  a  landing  time  to  each 
inbound  aircraft.  Once  the  estimated  landing  time  is  known,  the  speed  the  aircraft 
must  fly  needs  to  be  calculated.  This  cannot  be  done  directly  and  interpolation  using 
a  suitable  polynomial  approximation  is  used. 

This  memorandum  investigates  four  polynomials  and  examines  their  effectiveness 
at  providing  a  good  estimate  with  minimum  computation. 
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1  Introduction 


Air  traffic  movements  in  the  London  area  are  at  a  high  level  and  are  expected  to  increase 
Techniques  involving  the  use  of  computer  assistance  are  being  investigated  to  enable  air 
traffic  management  and  control  systems  to  handle  this  increase.  One  technique  in  par¬ 
ticular  being  studied  by  the  Terminal  Control  Systems  Development  Group  (TCSDG)  in 
AD4  Division  of  RSRE,  is  the  early  adjustment  of  the  speed  of  arriving  aircraft  so  that 
the  rate  of  flow  near  to  the  airports  is  closely  matched  to  landing  capacity  [lj. 

The  TCSDG  is  developing  a  Speed  Control  Advisor  (SCA)  [2,3] .  This  uses  an  algorithm 
which  includes  the  allocation  of  a  landing  slot  to  each  inbound  aircraft  based  on  a  specified 
runway  separation.  The  calibrated  airspeed  (CAS)  required  for  descent  in  order  that  the 
slot  time  is  achieved  needs  to  be  computed.  That  is,  given  the  time  over  which  an  aircraft 
must  fly  from  A  to  B  (usually  descent)  at  what  speed  must  it  travel? 

The  speed  value  needed  to  cause  an  aircraft  to  take  a  specified  time  to  fly  from  some 
point  A  to  another  point  B  along  a  specified  route  through  a  specified  wind  field  cannot 
be  calculated  directly.  The  reasons  for  this  are  that  aircraft  normally  operate  at  constant 
CAS  or  constant  Mach  Number  rather  than  at  constant  true  air  speed  (TAS),  that  different 
aircraft  types  perform  differently  especially  in  descent  and  turns,  and  that  wind  speed  and 
direction  can  vary  rapidly  with  change  of  altitude.  However,  the  time  taken  to  fly  from  A 
to  B  at  a  specified  CAS  can  be  calculated,  and  by  doing  this  calculation  for  several  speed 
values  and  using  some  form  of  interpolation,  the  required  speed  can  be  found. 

The  time  to  fly  from  A  to  B  can  be  expressed  as 


T  = 


i: 


where 

V7  =  V7,  +  w 

Vj  is  the  ground  velocity  vector,  V,  is  the  true  velocity  vector  and  W  is  the  wind  vector. 
V(  is  assumed  to  head  in  a  direction  such  that  Vg  points  along  the  route.  V,  is  a  function 
of  CAS  and  altitude,  or  of  Mach  number  and  altitude,  depending  on  the  phase  of  flight. 

A  set  of  procedures  for  performing  these  time  calculations  was  written  for  use  in  the 
TCSDG  real-time  system,  and  was  available  for  use  in  this  investigation. 

For  each  aircraft,  the  SCA  finds  the  time  taken  to  descend  at  maximum  speed,  minimum 
speed  and  preferred  speed  (the  speed  at  which  the  aircraft  would  like  to  descend)  These 
known  points  are  substituted  into  a  polynomial  approximation  to  give  an  estimate  (chosen 
speed)  for  the  required  CAS  (target  speed).  The  time  to  fly  the  route  portion  at  the  chosen 
speed  is  calculated  using  the  TCSDG  procedures,  and  is  compared  to  the  target  time.  If 
the  difference  is  not  within  acceptable  limits  then  further  iterations  are  needed.  As  each 
calculation  of  time  is  lengthy  it  is  important  to  find  a  polynomial  that  gives  a  good  estimate 
quickly. 

This  note  investigates  four  approximations  that  could  be  used  and  examines  their  effec¬ 
tiveness  at  providing  a  good  estimate  with  minimum  computation. 
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The  Approximations 


For  an  aircraft  flying  a  typical  descent  profile,  the  relationship  between  speed  and  time  is 
represented  by  the  actual  data  curve  in  Figs(l)-(4),(10)-(12).  (See  section  4  for  a  detailed 
description  of  descent  profile  used). 

Four  approximations  to  this  curve  have  been  tried.  In  the  derivations  that  follow  the 
notation  used  is 


U 

I. 

Ip 

Ic 

It 


maxcas  T j 

mincas  T, 

prefects  Tp 

chosen  speed  Tc 
target  speed  Tt 


time  to  fly  descent  at  maxcas 
time  to  fly  descent  at  mincas 
time  to  fly  descent  at  prefcas 
time  to  fly  descent  at  chosen  speed 
target  time 


where  I,  <  If  <  1/  and  Tf  <  Tp  <  T,  (since  it  takes  longer  to  fly  at  a  slower  speed) 


2.1  Approximation  I 

Initially,  only  two  known  points  were  used  ie  (If,Tf)(Is,Tr)  in  finding  a  suitable  approxi¬ 
mation.  The  first  assumption  made  was  that  the  relationship  between  speed  and  time  was 
inversely  proportional 

I  =  k/t  ( k  const) 

We  can  find  some  constant  k  such  that 

it -i,  =  k(i/Tf  -  i/r,) 

=  k(T ,  -  T,)!(T,Tf) 

=*  k  =  (If  -  I.)T,T.f(T,  *  Tf)  (1) 

Also,  given  that  Ic  is  proportional  to  Tc  then 

Ic  -  If  =h(l/T,-  l/7» 

=  k(Tf  -  Tt)/(TfT<) 

~IC  =  I,  +  k(T,  -  T,)/(TfT,) 

Substituting  for  k  from  (1)  gives 

/,  =  //  +  (//  -  It)(Tf  ~  Tt)TfT,/TfTt{T,  -  T<) 

=  It  +  (//  -  h)(T,  -  Tt)T,/Tt(T,  -  Tf) 

Fig(l)  shows  the  curve  this  equation  represents  plotted  against  the  data  being  modelled 

The  assumption  that  /  =  k/t  appears  to  provide  a  reasonable  model  and  this  approxima¬ 
tion  was  tested  along  with  the  others  to  see  what  sort  of  estimate  would  be  produced. 

The  use  of  only  two  points  is  not  ideal.  It  is  more  likely  that  a  better  fit  will  be  obtained 
with  more  points  and  as  three  are  being  calculated  already  (ie  preferred  speed  and  time), 
three  should  be  used 
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2.2  Approximation  II 


Approximation  II  is  an  extension  of  approximation  I  which  makes  use  of  the  three  points 
Speed  is  still  assumed  to  be  inversely  proportional  to  time  but  as  well  as  the  curve  shaping 
constant  k ,  two  further  constants,  a  and  b,  are  included  which  allow  the  curve  to  be  moved 
to  fit  the  data. 

0  +  b) 

We  have  three  points  ( Ij,Tf)(Ip,Tp )  and  (7,,  77)  which  when  substituted  into  (2)  give  the 
equations 


If  +  a  =  k/(Tf  +  b) 

Ip  +  a  =  k/{Tp  +  b) 

7,  +  a  =  k/{T.  +  6) 

These  equations  can  be  solved  to  give  the  following  values  for  6  ,  a  and  k. 

b  =  -  Tp)  +  IpTp{T,  -  77)  +  I,T,(TP  -  Tj) 

If(Tr  -  T,)  +  Ip(T.  -  77)  +  /,(7>  -  TP) 

_  Ip(b  +  Tp)~  I, (b  +  Tf) 

( Tf  -  T,) 

k  =  (7/  +  a)(T ;  4  6) 

which  are  substituted  into  (2)  to  give  an  approximation  of  speed  given  time.  The  three 
points  (If,Tf)(Ip,Tp)  and  (7,,T,)  must  be  distinct  in  order  to  provide  three  solvable 
simultaneous  equations. 

Fig(2)  shows  the  curve  this  equation  represents  plotted  against  the  data  being  modelled 
2.3  Approximation  III 

A  possible  numerical  analysis  solution  to  the  problem  of  fitting  a  function  to  some  data 
points  is  to  use  the  Lagrange  polynomial  (4,5] .  This  is  a  polynomial  that  can  be  constructed 
in  such  a  way  that  it  passes  through  each  known  data  point  and  is  of  degree  less  than  or 
equal  to  the  number  of  points  provided  minus  one.  For  example  if  we  know  three  points, 
the  Lagrange  polynomial  which  passes  through  these  points  will  be  quadratic  or  linear. 

2.3.1  Construction  of  tfae  Lagrange  polynomial 

Assume  that,  given  rt-f  1  distinct  points  t,  (t  =  0, 1, .  .  . ,  n)  on  the  interval  with  function 
values  /(<,)  =  /,  (t  =  0,1,..., n)  we  can  construct  a  polynomial  of  degree  <  n  which 
passes  through  the  n  +  1  points. 

The  Lagrange  method  is  to  begin  by  expressing  this  polynomial  as 

7(0  =  Lo(0/o  +  7-i(()/i  +  . . .  +  Ln(t)ft 

where  each  Ljt(()  (k  =  0, 1, . . . ,  n)  is  a  polynomial  of  degree  not  exceeding  n.  In  order  for 
the  polynomial  to  pass  through  the  n  +  1  points  given,  the  following  condition  must  hold 

7(0=  /.  (<  =  0,1,...,  n) 
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ie.  at  any  of  the  given  points,  the  Lagrange  polynomial  value  is  equal  to  the  actual  function 
value.  This  condition  holds  if 


* t .  |  1  *  =  «' 

}  0  k±i 


So  for  example 

I(t\)  —  Lo{tl)fo  +  L\(ti)fi  +  . . .  +  Ln(ti)Jn 
=  0*/o+l*/l+0*/2  +  ...+  0*/n 

=  /t 

The  family  of  polynomials  which  satisfy  the  requirements  above  are  of  the  form 

L  it)  =  (*  ~  <<>)(*  -<!)■■(<-  <*-»)(<  ~  <*+»)•  _lLz  t„) 

‘  (<*  -  t0)(tk  -  tk+1)  t„) 

Note  that  the  numerator  of  L*(t)  is  a  product  of  all  factors  of  the  form  (t  -  f, )  except 
(t  -  tfc).  Since  each  Lk(t)  is  a  polynomial  of  degree  n  and  since 

/(<)  =  io(‘)/o  +  il(‘)/l  +  ---+M0/« 

then  I(t)  will  be  a  polynomial  of  degree  <  n. 

2.3.2  Use  of  the  Lagrange  polynomial 

We  have  three  known  points  {1  f,Tf){IF,Tp)(I  s,Ta)  so  n  =  2  and  the  Lagrange  polynomial 
can  be  constructed  as  follows 

Ht)  =  L0[t)If  +  L1(t)Ip+Ln(t)I. 

where 

~  (r/-T,,.)(r/-rT5  -  (r,,-r,)(r,.-r,) 

,  _  (IrML- rp) 


i*(0  = 


(T,  -  r/)(T,  -  Tp) 


Fig(3)  shows  the  curve  this  equation  represents  plotted  against  the  data  being  modelled. 
It  is  clearly  a  reasonable  fit  and  should  provide  a  good  estimate  of  the  speed  given  time. 


2.4  Approximation  IV 

The  fourth  approximation  assumes  that  the  speed/time  relationship  can  be  represented 
by  an  inverse  cubic. 

'(0  =  7  +  T?  +  ji  (3) 


We  have  three  points  (If,  Tj)(Ir,  Tp)  and  ( I,,T ,)  which  when  substituted  into  the  inverse 
cubic  give  the  equations 


1/  =  ao/Tf  +  ai/Tf  +  at/Tj 

(4) 

/p  =  aQ/Tp  +  ai/Tp  +  at/Tp 

(5) 

I,  =  ao/T,  +  ai/T?  +  at/T, 

(6) 

These  can  be  solved  to  give  values  for  ao  ,  ai  and  at  which  are  substituted  into  (3)  to 
give  an  approximation  for  speed  given  time.  The  three  points  ( If,T/)(Ip,Tp )  and  ( I,,T„ ) 
must  be  distinct  in  order  to  provide  three  solvable  simultaneous  equations. 
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2.4.1  Calculation  of  ao  ,  aj  ,  a2 


We  can  write  eqns  (4)-(6)  as 

(  If  \  (  1/77  i/Tj  l/T } 

/,  =  l/T,  1  /Tj  1/Tj 

\  I,  )  \  l/T,  l/T j  l/T f 

represented  by 

7=  Ta 

and  so  if  T1  is  the  inverse  of  T  then 

T'J  =  T'Ta  =  a 


1  ao 

t  °2 


The  inverse  of  T  is  defined  as 

T'  =  adj  T/  det  T 

where  adj  T  is  the  transpose  of  the  matrix  of  cofactors  of  T  and  det  T  is  the  determinant 
of  T.  As  long  as  T/  /  /  T,  then  det  T  will  not  equal  zero.  From  this  we  can  find 

ao  ,  Gi  and  a2 

Let  B  be  the  matrix  of  cofactors  of  T.  Then 


B  = 


f JT}  TjT? 


T.  T3  +  T$T.  T.Va 


Tp  +  Tfr; 


T,.T a 
1 

t,t; 

1 


'  T/Tj 


3  + 


1 

r, t? 

_ )_ 

TjT} 

TjT?  ’ 


TjT. 


t;t»t; 


TjT, 

tJt,  TjT}  -  Tfr„ 

(  Tj(Tp  -  T,)  T/(T;  -  T*)  T/(T.T(?  -  T}TP) 

Tj(T,  -  T/)  Tj(Tj  -  Tj)  Tj(TfTj  -  TjT.) 

V  T*{Tf  -  T,)  Tj(Tj  -  Tj)  Tj(TPTj  -  TjT;) 


The  adj  T  is  the  transpose  of  the  matrix  of  cofactors,  that  is  the  transpose  of  B. 
Hence  adj  T  — 

T}{T,-T,)  Tj(T,  -  Tj)  Tj(T,-Tp)  \ 

Tj(Tj  -  Tj)  Tj(Tj-Tj)  Tj{Tj  -  Tj) 

Tj(T,Tj  -  TjTp)  Tj(TfTj  -  TjT,)  Tj(TpTj  -  TjTj)  J 
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TjTjTj 

The  inverse  of  T  is  equal  to  adj  T j det  T 
det  T  = 


TjTjTj 


We  have 
and  so 


a0  = 


Tj(Tp  -  T,)  +  T2(T.  -  Tf )  +  T‘(Tf  -  Tp) 
T'l  =  a 

I/Tf{TP  -  T ,)  +  IpTj(T,  -  Tf)  +  I,Tj(T,  -  Tp) 
Tj(Tp  -  T,)  +  Tj(T,  -  Tf)  +  Tj(T,  -  Tp) 


<J2  = 


IfTj(Tj  -  Tj)  +  IpTj(Tj  -  Tj)  +  I,Tj(Tj  -  Tj) 

Ql  T}(TP  -  T,)  +  Tj(T,  -  Tf)  +  Tj[T,  -  T„) 

_  IfT}{T,T*  -  TjT,)  +  IPTj(T,Tj  -  TjT,)  +  I,Tj(TpTj  -  Tpj) 


T}(T,  -  T,)  +  Tj ( T,  -  T,)  +  T}(T{  -  Tp) 


Fig(4)  shows  the  curve  represented  by  the  inverse  cubic  compared  to  the  actual  data 
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3  Running  the  experiment 


The  aim  of  this  investigation  was  to  find  the  best  approximation  modelling  the  relationship 
between  speed  and  time.  It  is  possible  to  calculate  the  time  taken  to  fly  a  descent  given 
the  speed  by  using  a  series  of  complex  prediction  procedures.  This  process  is  lengthy  so 
the  best  approximation  will  be  one  that  takes  minimum  calls  of  these  procedures  to  deliver 
a  reasonable  estimate  (ie.  one  within  specified  limits)  of  speed  given  time  . 

One  program  for  each  algorithm  was  written  in  CORAL  66  The  programs  were  basically 
the  same  but  had  slight  variations  according  to  the  algorithm  being  used.  The  programs 
included  calls  to  the  prediction  procedures  which  calculate  the  time  to  fly  the  descent 
given  the  speed.  The  route  portion  chosen  was  of  an  aircraft  approaching  Heathrow  from 
the  North  descending  from  a  height  of  31,000ft  to  8,000ft  (stack  base  level)  The  descent 
profile  consisted  of  a  cruise  section  at  the  Mach  equivalent  of  290k  followed  by  a  descent 
at  the  chosen  speed.  This  descent  was  subject  to  minimum  calibrated  airspeeds  at  defined 
levels  illustrated  in  Fig(5). 

Initially,  the  times  taken  to  fly  the  route  portion  at  max  speed,  min  speed  and  preferred 
speed  were  found,  providing  the  three  points  (360, 1 153), (210, 1569),  (260,1333)  These 
were  used  to  find  the  constants  in  each  algorithm  A  target  time,  Tt  .  was  chosen  ie 
the  time  the  aircraft  had  to  fly  the  descent  and  an  estimate  for  the  speed  required  to  fly 
the  descent  was  found  by  substituting  the  target  time  into  each  approximation  Having 
obtained  this  value  for  the  speed,  the  investigation  proceeded  to  see  how  reasonable  it 
was. 

The  prediction  procedures  were  called  to  calculate  the  actual  time  an  aircraft  would  take 
given  this  estimated  speed  If  the  difference  between  actual  time  and  target  time  was  less 
than  5  secs,  the  estimated  speed  was  returned  as  the  chosen  speed,  Ic  ,(the  s[  eed  needed  to 
fly  the  descent).  If  the  difference  was  greater  than  5  secs,  further  iterations  were  required 

3.1  Further  iterations 

The  calculation  of  the  time  to  fly  the  route  portion  at  estimated  speed  provided  an  extra 
point  (Ic,Te).  Since  Tc  was  nearer  to  the  target  time  than  at  least  one  of  the  end  points, 
it  was  used  in  the  further  iterations  to  obtain  a  better  estimated  speed 

3.1.1  Approximation  I 

The  obvious  method  to  use  with  this  type  of  approximation  was  to  shorten  the  interval 
over  which  the  interpolation  was  taking  place  This  was  achieved  using  the  secant  method 


3.1.2  Approximation  II 

The  method  used  with  this  approximation  shortened  the  interval  over  which  interpolation 
was  taking  place.  In  this  case  the  interval  was  over  I,  to  Ij  with  Ip  in  between.  Moving 
of  the  end  points  had  to  take  into  account  all  three.  The  extra  point  ( /C,TC )  was  taken  to 
be  the  mid  point  of  the  next  three  points  used  since  a  better  fit  was  achieved  The  end 
points  were  then  picked  accordingly. 
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3.1.3  Approximation  III 


The  extra  point  (IC,TC)  was  used  in  addition  to  the  three  original  points,  to  find  the  next 
order  Lagrange  approximation  eg.  one  with  degree  <  4.  This  gave  a  function  that  passed 
through  the  four  known  points. 

3.1.4  Approximation  IV 

The  method  used  with  this  inverse  cubic  approximation  was  the  same  as  for  approximation 

II 


3.2  Output 

Each  program  was  run  with  target  times  ranging  from  1165  to  1548  secs  which  produced 
speeds  between  214  and  329  knots. 

Target 

time  1165  1195  1236  1278  1319  1361  1403  1444  1486  1527  1548 

Speed  329  309  290  276  264  253  244  236  228  219  214 

The  target  times  were  chosen  to  give  a  fairly  even  spread  over  the  range  between  Tj  and 
T. 

The  output  recorded  for  each  target  time  was 

•  estimated  speed 

•  actual  time  for  descent  at  estimated  speed 

•  final  chosen  speed  and  time 

•  number  of  iterations  required 

The  algorithms  were  compared  by  looking  at  the  number  of  iterations  required  to  achieve 
a  good  estimate  This  shows  how  often  the  prediction  procedures  were  called  giving  an 
indication  of  how  good  the  algorithm  was  and  how  closely  it  modelled  the  speed  time 
relationship 

The  Extended  inverse  proportional,  the  Lagrange  and  the  Inverse  cubic  approximations 
were  further  examined  to  see  how  the  choice  of  the  third  point  (initially  the  preferred  cas) 
affected  the  reasonableness  of  the  estimate 


4  Results 

A  “reasonable  estimate”  of  the  speed  was  defined  as  one  whose  associated  time  to  fly  the 
route  portion  was  within  5secs  of  the  target  time  Every  computation  of  an  estimated 
speed,  including  the  first,  was  considered  to  be  an  iteration.  The  number  of  iterations 
required  for  each  algorithm  are  shown  in  Fig(6) 

From  this  it  can  be  seen  that  the  Extended  inverse  proportional,  Inverse  cubic  and  La¬ 
grange  algorithms  provide  good  approximations,  taking  at  most  2  iterations  to  give  a 


7 


reasonable  estimate  The  Inverse  proportional  algorithm  provides  a  fair  approximation, 
taking  3  iterations  at  most  speeds  but  less  near  each  end  point. 

It  can  be  seen  that  the  algorithms  using  three  points  provide  a  better  model  than  the 
Inverse  proportional,  especially  around  the  preferred  speed  mark.  This  is  to  be  expected 
since  these  three  polynomial  approximations  are  constructed  in  such  a  way  that  they  pass 
through  the  three  original  points  (//,  Tj)(Ir,  Tp)(If,  T,).  The  implication  of  this  is  that  the 
choice  of  preferred  speed  has  an  effect  on  the  accuracy  of  the  model,  the  question  being 
“How  much  ?”. 

In  order  to  study  this  effect,  the  experiment  was  run  with  preferred  speeds  of  330k,  275k 
and  320k.  Figs  (7), (8), (9)  show  the  difference  between  the  number  of  iterations  required 
for  each  algorithm. 

The  fig_.es  show  that  for  speeds  around  the  preferred  speed  and  two  end  points,  the 
estimate  is  close  to  the  target.  Between  these  points,  the  estimates  get  worse  as  the 
models  become  less  accurate. 

With  a  preferred  speed  of  230k  ( Fig ( 7) ) ,  the  algorithms  are  providing  reasonable  models, 
taking  mostly  two  iterations  although  the  Inverse  cubic  and  Lagrange  take  three  iterations 
at  speeds  just  below  340k. 

With  a  preferred  speed  of  275k  (Fig(8)),  the  actual  midway  point  between  the  end  speeds, 
all  three  algorithms  provide  good  approximations  taking  no  more  than  two  iterations. 

A  preferred  speed  of  320k  (Fig(9) )  gives  rise  to  the  least  accurate  estimates  The  Lagrange 
algorithm  is  particularly  affected  by  this  high  preferred  speed  taking  three  or  four  iterations 
to  converge  on  speeds  away  from  the  end  points.  The  Inverse  cubic  and  Extended  inverse 
proportional  are  also  affected  but  not  to  such  a  great  extent 

These  results  are  illustrated  by  Figs(10),(l  1)  and  (12)  which  show  the  graphs  of  actual  data 
along  with  the  graphs  of  the  Extended  inverse  proportional.  Inverse  cubic  and  Lagrange 
algorithms  with  differing  preferred  speeds.  These  demonstrate  the  effectiveness  of  each 
algorithm  to  provide  a  model  of  the  actual  data. 


5  Conclusions 

The  four  possible  relationships  considered  between  speed  and  time  were 

1.  inverse  proportional  ie  I  =  k/t 

2.  extended  inverse  proportional  ie  I  +  a  =  k/(t  +  ft) 

3.  quadratic  (modelled  with  three  points)  using  the  Lagrange  polynomial 

4.  inverse  cubic  ie  I  =  ao/t  +  ai/tJ  +  aj/t3 

The  algorithms  were  set  up  using  three  points  of  data,  the  minimum,  maximum  and 
preferred  speeds  and  their  associated  time  of  descents.  A  target  time  was  input  into  each 
algorithm  and  an  estimate  for  speed  found.  The  time  taken  to  descend  at  this  estimated 
speed  was  then  compared  with  the  target  time.  A  reasonable  estimate  was  defined  to  be 
one  within  5secs  of  the  target  time. 

The  inverse  proportional  relationship  provided  a  rough  model  and  gave  a  reasonable  esti¬ 
mate  within  three  iterations. 
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The  Lagrange  polynomial  provided  a  good  model,  giving  a  reasonable  estimate  within 
two  or  three  iterations  using  the  low  to  centre  preferred  speeds  However,  the  choice  of  a 
high  preferred  speed  had  a  large  effect  on  the  first  estimate  and  up  to  four  iterations  were 
required. 

The  inverse  cubic  also  provided  a  good  model,  giving  a  reasonable  estimate  within  two  or 
three  iterations  using  any  of  the  preferred  speeds. 

The  extended  inverse  proportional  provided  the  best  model  of  the  relationship  between 
speed  and  time.  The  algorithm  took  the  same  number  of  iterations  or  less  than  the  other 
approximations  and  was  not  affected  unduly  by  the  change  in  preferred  speed.  The  fit 
of  the  extended  inverse  proportional  to  the  actual  data  points  was  the  closest  of  all  the 
algorithms  (Fig(12)),  and  had  the  computational  advantage  that  its  constants  were  smaller 
than  those  for  the  inverse  cubic. 

This  algorithm  has  now  been  incorporated  into  the  Speed  Control  Adviser.  It  has  been 
found  that  the  target  CAS  is  normally  given  directly  (to  within  required  tolerance  lev¬ 
els)  although  occasionally  two  or  three  iterations  are  necessary  The  extended  inverse 
proportional  algorithm  is  therefore  working  effectively  and  keeping  processing  time  to  a 
minimum. 
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Figure  8:  Preferred  CAS  =  275K 
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Figure  9:  Preferred  CAS  =  320K 


P  =  Inversely  proportional  L  =  Lagrange 
E  =  Extended  1/T  C  =  Inverse  cubic 
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